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ABSTRACT 

We present the results of a study on the kiloparsec-scale radio emission in the 
complete flux density limited MOJAVE sample, comprising 135 radio-loud AGNs. 
New 1.4 GHz VLA radio images of six quasars, and previously unpublished im- 
ages of 21 blazars, are presented, along with an analysis of the high resolution 
(VLA A-array) 1.4 GHz emission for the entire sample. While extended emission 
is detected in the majority of the sources, about 7% of the sources exhibit only 
radio core emission. We expect more sensitive radio observations, however, to 
detect faint emission in these sources, as we have detected in the erstwhile "core- 
only" source, 1548+056. The kiloparsec-scale radio morphology varies widely 
across the sample. Many BL Lacs exhibit extended radio power and kiloparsec- 
scale morphology typical of powerful FRII jets, while a substantial number of 
quasars possess radio powers intermediate between FRIs and FRIIs. This poses 
challenges to the simple radio-loud unified scheme, which links BL Lacs to FRIs 
and quasars to FRIIs. We find a significant correlation between extended radio 
emission and parsec-scale jet speeds: the more radio powerful sources possess 
faster jets. This indicates that the 1.4 GHz (or low frequency) radio emission 
is indeed related to jet kinetic power. Various properties such as extended ra- 
dio power and apparent parsec-scale jet speeds vary smoothly between different 
blazar subclasses, suggesting that, at least in terms of radio jet properties, the 
distinction between quasars and BL Lac objects, at an emission-line equivalent 
width of 5A is essentially an arbitrary one. While the two blazar subclasses 
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display a smooth continuation in properties, they often reveal differences in the 
correlation test results, when considered separately. This can be understood if, 
unlike quasars, BL Lacs do not constitute a homogeneous population, but rather 
include both FRI and FRII radio galaxies for their parent population. It could 
also just be due to small number statistics. We find that the ratio of the radio 
core luminosity to the /c-corrected optical luminosity (R v ) appears to be a better 
indicator of orientation for this blazar sample, than the traditionally used radio 
core prominence parameter (R c ). Based on the assumption that the extended ra- 
dio luminosity is affected by the kiloparsec-scale environment, we define the ratio 
of extended radio power to absolute optical magnitude (L ext /M a b s ) as a proxy for 
environmental effects. Trends with this parameter suggest that the parsec-scale 
jet speeds and the parsec-to-kiloparsec jet misalignments are not affected by the 
large-scale environment, but are more likely to depend upon factors intrinsic to 
the AGN, or its local parsec-scale environment. The jet speeds could, for in- 
stance, be related to the black hole spins, while jet misalignments could arise due 
to the presence of binary black holes, or kicks imparted to black holes via black 
hole mergers, consistent both with radio morphologies resembling precessing jet 
models observed in some MOJAVE blazars, and the signature of a 90° bump 
in the jet misalignment distribution, attributed to low-pitch helical parsec-scale 
jets in the literature. We suggest that some of the extremely misaligned MO- 
JAVE blazar jets could be "hybrid" morphology sources, with an FRI jet on 
one side and an FRII jet on the other. Finally, it is tempting to speculate that 
environmental radio boosting (as proposed for Cygnus A) could be responsible 
for blurring the Fanaroff-Riley dividing line in the MOJAVE blazars, provided a 
substantial fraction of them reside in dense (cluster) environments. 

Subject headings: galaxies: active — quasars: general — BL Lacertae objects: 
general — radio continuum: galaxies 



1. INTRODUCTION 

Radio-loud active galactic nuclei (AGNs) that are characterised by extreme variability 
in their radio cores, high and variable polarization, superluminal jet speeds and compact 
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radio emission, are collectively referred to as 'blazars' (lAngel fc Stockman Il980h . Blazars 
comprise flat-spectrum radio-loud quasars and BL Lac objects. Their extreme characteristics 
have commonly been understood to be a c onsequence of relatiyistic beaming effects due to 
a fast jet aligned close to our line of sig ht falandford fc Konigll[l979l ). 



The radio-loud unified scheme postulates t hat quasars and BL L acs are the beamed 
end-on counterparts of the Fanaroff-Riley (FR. iFanaroff fc Riley 1974 ) type II and type I 
radio galaxies, respectively ( lUrry fc Padovanilll995l ). The relatively lower radio luminosity 
FRI radio galaxies are referred to as "edge-darkened" since their brightest lobe emission 
lies closer to the radio cores and fades further out, while the higher radio luminosity FRII 
radio galaxies are "edge-brightened" because of the presence of bright and compact hot 
spots where the kpc-scale jets terminate. This definition turns out to be useful for FRI-FRII 
classification when hot spots are not clearly delineated (cf. §3.2). 

The Fanaroff-Riley dichotomy has been proposed to arise due to differences in one 



or more of the following pro perties: host galaxy environmen t (IPrestage fc Peacock! 11988 



and jet-me dium interaction (IBickne 



mode ( Baum et al. 



1995 



hole m ass fjGhisellini fc Celottil 120011), black hole sp in ( iMeierl Il999l ) . accretion rate and/or 



ll Il995l ). jet compo sition ([Reynolds et al.lll996l ). black 



Ghisellini &: Celotti 



2001 



Many recent findings are, however, 
posing c hallenges for the sta ndard unified sc heme. These incl ude the discovery of FRI 
quasars (IHeywood et al.l 120071 ). FRII BL Lacs (ILandt et al.ll2006l ). and "hybrid" radio mor- 
phology sources with an FRI jet on one side of the core and an FRII jet on the other 
(IGopal-Krishna fc Wiital 120001 ) . One of the primary cited differences between quasars and 
BL Lacs has been the presence of s trong, broad e missi on lin es in the quasar sp ectra and 
their apparent absence in BL Lacs. IStickel et al.l ( 1l99ll ) and IStocke et al.l (119911 ) have de- 
fined the distinction between BL Lacs and quasars at an e mission line equivalent width 
of 5 A However, this d istinction has been questioned (e.g., IScarpa fc Falomol 119971 ; lUrry 



1999 



Landt 



Miller et al 



et al. 



19781 : 



200 4J), and it is known that some BL Lacs have a broad-line r egion (e.g. 



Stickel et al.1ll99ll : IVermeulen et al.lll995l : ICorbett et al"1ll996h 



In this paper, we examine the unified scheme and FR dichotomy in the MOJAVtlj 
sample of blazars. MOJAVE is a long-term program to monitor radio brightness and po- 
larization variations in the jets associated with active galaxies vi sible in the northe rn sky 
on parsec-scales with Very Long Baseline Interferometry (VLBI) (ILister et al.ll2009al ). The 
MOJAVE sample consists of 135 sources satisfying the following criteria: (1) J2000.0 dec- 
lination > —20°; (2) galactic latitude |6| > 2.5°; (3) VLBA 2 cm correlated flux density 



2 Monitoring Of Jets in Active galactic nuclei with VLBA Experiments. 
http://www.physics.purdue.edu/MOJAVE/ 
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exceeding 1.5 Jy (2 Jy for declination south of 0°) at any epoch between 1994.0 and 2004.0. 
Of the 135 AGNs, 101 are classified as quasars, 22 as BL Lac objects, and eight as radio 
galaxies of mostly the FRII-type. Four radio sources are yet to be identified on the basis 
of their emission line spectra and have no redshift information. Overall, eight sources (four 
BL Lacs and four unidentified) have unknown or unreliable redshifts. Based on the syn- 
chrotron peak in the spectral energy distributions (SEDs), BL Lacs h ave been divided into 
low, h igh, and intermediate energy pe aked classes (LBLs, HBLs, IBLs, iPadovani fc Giommi 
199.4 lLaurent-Muehleisen et al.lll999l ). All but three MOJAVE BL Lacs are classified as 

The three BL Lacs, viz., 0422+004, 1538+149 and 



Nieppola et al 



200 



J NED). 



LBLs (see 

1807+698, are classified as IBLs. There are no HBLs in the MOJAVE sample. 

The compact flux density selection criteria of the MOJAVE sample biases it heavily to- 
ward highly beamed blazars with high Lorentz factors and small viewing angles. While 
relativistic boosting effects are likely to dominate the source characteristics, the exten- 
sive multi-epoch and multi-wavelength data available for the MOJAVE sample on both 
parsec- and kiloparsec-scales, provides us with unique constraints to test the unified scheme. 
Throughout the paper, we adopt the cosmology in which H =71 km s -1 Mpc -1 , f2 m =0.27 
and f2A=0.73. Spectral index, a, is defined such that, flux density S u at frequency u, is 
S u oc v~ a . 



2. DATA REDUCTION AND ANALYSIS 

We observed seven MOJAVE quasars at 1.46 GHz with the Very Large Array (VLA, 
Napier et al.lll983h in the A-array configuration (typical synthesized beam ~ 1.5") on June 



30, 2007 (Program ID:AC874). Si xty MOJAVE s ource s were previously observed by us with 



the VLA A-array and presented in lCooper et al.l (120071 ). Data for the remaining sources were 



reduced directly using archival VLA A-array data, or obtained through published papers. 
The data reduction was carried out following standard calibration and reduction procedures 
in the Astronomical Image Processing System (AIPS). 

Our new observations included some new expanded VLA (EVLA) antennas and the 
observations were made in a spectral line (4-channel pseudo-continuum) mode, which resulted 
in a total effective bandwidth of (2 x 21.875 =) 43.75 MHz. Four m 9 minute scans of each 
source were interspersed with 1 minute scans of a suitable phase calibrator. 0713+438 was 
used as the bandpass calibrator for the experiment. After the initial amplitude an d phase 



calib ration, AIPS tasks CALIB and IMAGR were used iteratively to self-calibrate (ISchwab 



1980) and image the sources. The resultant rms noise in the maps was typically of the order 



of 0.06 mJy beam 1 . The quasar 1124—186 was observed with an incorrect source position, 
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Fig. 1.— 1.4 GHz VLA images of the quasars 0730+504 (Left) and 0805-077 (Right). The 
contours are in percentage of the peak surface brightness and increase in steps of 2. The lowest 
contour levels and peak surface brightness are (Left) ±0.042, 672 mJy beam -1 and (Right) ± 0.010, 
1.43 Jy beam -1 . 

which resulted in significant beam smearing. We instead obtained archival VLA data for 
this source. The new radio images of the six sources are presented in Figures [H [21 and El 

Since most of the MOJAVE sources have highly compact radio structures, most of the 
archival data comprised of snapshot observations where the sample sources were observed 
as phase calibrators. By choosing archival datasets with exposure times >10 minutes, we 
were able to obtain maps with typical rms noise levels of ~ 0.15 mJy beam -1 . We present 
previously unpublished images of 21 blazars in Appendix A. 

A compilation of the basic parameters for each source is given in Table [TJ The integrated 
flux densities for the cores were obtained in AIPS using the Gaussian-fitting task JMFIT, 
while the total flux densities were obtained by putting a box around the source, using the 
AIPS verbs TVWINDOW and IMS TAT. The extended radio flux density was obtained by 
subtracting the core flux density from the total radio flux density. 

All the maps were created with uniform weighting using a ROBUST parameter of in 
the AIPS task IMAGR. We also obtained core and lobe flux densities from maps with different 
weighting schemes (by using ROBUST parameters —5 and +5, respectively) for a fraction 
of the sample. We found that the integrated core flux density estimates differed typically by 
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Fig. 2 — 1.4 GHz VLA images of 1036+054 (Left) and 1045-188 (Right). The contours are in 
percentage of the peak surface brightness and increase in steps of 2. The lowest contour levels and 
peak surface brightness are (Left) ±0.021, 892 mJy beam -1 and (Right) ± 0.042, 724 mJy beam -1 . 



less than 1%, between different weighting schemes. The extended flux density values differed 
typically by < 2%, with a large majority of sources showing less than a 5% difference. Only 
in a handful of sources, where there seemed to be an unresolved compact component close to 
the core (e.g., 1038+064, 1417+385), was the difference between dif ferent weighting sche mes 
significant (10%— 20%). However, following our radio galaxy study (IKharb et al.ll2008bh . we 
have found that nearly 15%— 20% of the extended flux density could be getting lost in A-array 
observations, compared to combined-array observations. Therefore, we conclude that the lack 
of short spacings in the A-array observations is far more detrimental to the determination 
of accurate extended flux values than (not) adopting different weighting schemes on the 
A-array data to obtain core and extended flux densities. The weighting scheme approach 
does suggest that the errors in the extended flux density values are typically of the order 
of 2%— 5%, but could be of the order of 10%— 20% for sources where a compact component 
close to the core is not clearly resolved. 




12 15 47.6 47.4 47.2 47.0 46.8 46.6 46.4 46.2 46.0 12 22 23.4 23.2 23.0 22.8 22.6 22.4 22.2 22.0 21.8 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 



Fig. 3.— 1.4 GHz VLA images of 1213-172 (Left) and 1219+044 (Right). The contours are in 
percentage of the peak surface brightness and increase in steps of 2. The lowest contour levels and 
peak surface brightness are (Left) ±0.021, 1.66 Jy beam -1 and (Right) ± 0.042, 549 mJy beam -1 . 



3. RESULTS 



3.1. Extended Radio Power 



The 1.4 GHz extended radio luminosities for the MOJAVE sources are plotted against 
redshift in Figure HI Sources with no discernible extended emission are represented as upper 
limits. The solid lines indicate the FRI— FRII di vide (extrapolate d from 178 MHz to 1.4 
GHz a ssuming a spectral index, a = 0.8), following iLedlow fc Owenl ( 119961 ) and lLandt et al. 



(120061 ). The right hand panel of Figure H] demonstrates the close relation between the radio 
core and extended luminosity (Table Using the partial correlation regression analysis 
routine in IDL (P .CORRELATE), we found that the linear correlation between \ogL core 
and \ogL ext is strong even with the effects of luminosity distance (log-Dz,) removed (partial 
correlation coefficient, rxvz, = 0.303, t statistic = 3.41, two-tailed probability that the 
variables are not correlated^!, p = 0.0009). We note that the 18 BL Lacs (ones with redshift 
information) alone fail to show a correlation between logL core and logL ext , when the effects 
of luminosity distance are removed. The implication of this finding is discussed ahead in 



3 Calculated using the VassarStats statistical computation website, 

http://faculty.vassar.edu/lowry/VassarStats.html 
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§3.3. 
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Fig. 4. — (Left) 1.4 GHz extended luminosity versus redshift. (Right) 1.4 GHz core versus extended 
luminosity. Open and filled circles denote quasars and BL Lacs, respectively, while open squares 
denote radio galaxies. Core-only sources are represented as upper limits. 

The salient features of the extended radio emission in the MOJAVE blazars are: 



(i) Six of the 22 BL Lacs (~ 27%) have FRII radio powers {logL l J t > 26). Four of these 
(viz., 0235+164, 0716+714, 0808+019, 2131-021) also have hot spots like quasars. 

(ii) Five BL Lacs (~ 23%) fall in the FRI/II power range (24.5 < logL^ < 26). Two of 
these (viz., 0823+033, 1803+784) also appear to have hot spots. 

(iii) Seven BL Lacs have extended luminosity \ogL\^ t < 24.5 and could be regarded as 
the true beamed counterparts of FRI radio galaxies. However, three of these (viz., 
0754+100, 0851+202, 1807+698) also appear to exhibit hot spot-like features. Hot 
spots may also be present in all four BL Lacs with no redshift information. Overall, 
nearly 60% of the MOJAVE BL Lacs appear to have hot spots. The hot spots in some 
BL Lacs are however not as bright or compact as those observed in quasars. 

(iv) Excluding the upper limits, 22 quasars (~ 22%) fall in the FRI/II power range. 

(v) 10 sources (~ 7%, 9 quasars, 1 BL Lac) do not show any extended emission — these 
are discussed in §3.4. 
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Based on the extended emission at 1.4 GHz, we can conclude that a substantial frac- 
tion of the MOJAVE BL Lacs have both radio powers and radio morphologies like FRIIs 
or quasars. A substantial fraction of the MOJAVE quasars lie in the intermediate (FRI/II) 
luminosity range. These results are consistent with a number of p r eviou s radio studies. 
Using a large sample of radio core-dominated AGNs, iMurphv et al.l (119931 ) observed that 
many high redshift (z > 0.5) BL Lacs lie above the FRI /FRI I luminosity divi s ion. Using 
a sample of 17 BL Lacs, mostly belonging to the 1-Jy sample, iKollgaard et al.l (119921 ) sug- 
gested that about 12% — 30% of BL Lacs in a radio flux-limited sample may be bona fide 
FRIIs. Based on exten ded radio emission and strong emission lines at one or more epochs, 



Rector fe Stockd (120011) co ncluded that many radio-selected BL Lacs belonging to the 1 Jy 



sample (IStickel et al.lll99lf) cannot be beamed FRIs, but are more likely to be beamed FRIIs. 



Cara fe Lister! (120081 ) found similar intrinsic parent luminosity functions for the MOJAVE 



sample (consistent with FRIIs) irrespective of whether the BL Lacs were included or not. 



3.2. Radio Morphology 



We observe from Figures [TH3] that all but two sources show two-sided radio structures. 
The variety of r adio structures obs erved here are reasonably representative of the entire 
sample (e.g., see ICooper et al.l 120071'). a s well as those previously observed in other quasar 
surveys (e.g., 



Gower fe Hutchings 



1984 ; lAntonucci fe Ulvestadl Il985l : IPearson fe Readhead 



19881 : IMurphv et al.l Il993l ). The quasars often show two-sided radio structures with hot 
spots on one or both sides of the core. However, one-sided morphologies with no discernable 
compact or diffuse emission on the other side of the core, are also common. 

While determining if the blazars had an FRII type radio morphology, we relied initially 
on the presence of one or more compact hot spots at the leading edge of the radio lobes. 
When no clear hot spots were visible on either side of the cores, we resorted to the traditional 
"edge-brightened" definition for FRII sources, i.e., when the brightest radio emission was 
furthest from the core, we classified the source as an FRII. The remaining sources with no 
clear hot spots, and the brightest extended emission closest to the cores, were classified as 
FRI types. Even then, it was sometimes difficult to consign the sources to FRI or FRII 
classes, and we have listed more than one morphology for some sources in Table HJ 

Many of the MOJAVE sources show distinctly curved kiloparsec-scale jets. Straight jets 
are a rarity in the sample. Many MOJAVE quasars exhibit large parsec-to-kiloparsec jet 
misalignments. While any intrinsic curvature in the jets is likely to be highly exaggerated 
by projection effects in these low-viewing angle blazars, many sources seem intrinsically 
distorted. We discuss in §4.1 how the MOJAVE selection criteria might be preferentially 
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picking up bent jets. Apart from highly curved jets, many sources show distinct hot spot-like 
features both closer to the core and at the jet terminati on points, similar to the wide angle 
tail (WAT) radio galaxies (e.g., O'Donoghue et al"1ll99ol ). 



Four MOJAVE blazars have not yet been identified as quasars or BL Lacs (Tabled]). 
A new VLA image of one of the unidentified sources (viz., 1213—172) is presented in Fig- 
ure US Another unid entified source (0648—165) has a similar core-dominant radio morphology 
( iCooper et al.ll2007l ). A third source (2021+317) has a distinct core-halo morphology resem- 
bling a BL Lac object (see App endix A), while the fourth (0446+112) has a straight jet with 
a possible hot spot at the end ( ICooper et al.l 120071 ). 



3.3. Parsec-scale Jet Speeds and Extended Luminosity 
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Fig. 5. — (Left) 1.4 GHz core luminosity versus the apparent jet speed. The aspect curve assumes 
7 = 52, Li n t = 5 x 10 24 and p ~2. (Right) 1.4 GHz extended luminosity versus the apparent jet 
speed. Open and filled circles denote quasars and BL Lacs, respectively, while open squares denote 
radio galaxies. Core-only sources are represented as upper limits. 



The MOJAVE and 2 cm survey programs have provided apparent parsec-scale jet 
speeds ((3 app ) for nearly the entire MOJAVE sample, by monitoring structura l changes over 



timescales of approximately a decade. Thes e have b een previously tabulated by lKellermann et al. 
( 120041 ) . and most recently by iLister et al.l (l2009bl ). For our analysis below, we use a j3 app 
value that represents the fastest moving feature in that source. 



The parsec-scale apparent jet speeds and kiloparsec-scale radio core luminosity (Fig. E]) 
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are related by standard beaming relations. The aspect curve in Figure El assumes: L = 
Lint x 8 P , fiapp = sin 6/(1 — (3 cos 6), where the Doppler factor, 5 = 1/(7 (1 — (3 cos 6)) and 
= v/c. The best-fit values for the curve are 7 = 52, Li nt = 5 x 10 24 and p ~ 2. The relation 
between parsec-scale apparent jet spee ds and parsec-scale core l uminosity for the MO JAVE 



sample is discussed in greater detail bv lCohen et al.l ( 120071 ). and iLister et al.l ( 12009bl ). 



We find that the parsec-scale apparent jet speeds are correlated significantly with the 
extended radio luminosity (Fig. El Table [2]). Partial regression analysis shows that the linear 
correlation between \ogL ext and /3 app with the effects of luminosity distance removed, is still 
highly significant (txy.z = 0.273, t statistic = 3.05, two-tailed probability, p = 0.0028). 
This implies that more radio powerful sources have faster radio jets. Since extended radio 
luminosity is correlated with the core luminosity (Fig. HJ Table and the core luminosity 
is related to the apparent jet speed, we tested the correlation with a partial correlation 
test between extended radio luminosity and apparent jet speeds, after removing the effects 
of radio core luminosity. This too yielded a statistically significant correlation between 
extended radio luminosity and parsec-scale apparent jet speeds (r X y.z = 0.246, t statistic = 
2.72, p = 0.0075). 

The significant implication of this result is that faster jets are launched in AGNs with 
larger kiloparsec-scale lobe luminosities. It would therefore appear that the fate of the AGN 
is decided at its conception. This result undermines the role of the kiloparsec environment 
on the radio power of a given source. It also indicates that the 1.4 GHz extended emission 
is indeed related t o jet kinetic power, contrary to some suggestions in the literature (e.g., 



Birzan et al.ll2004l ). Given that there is a large overlap in radio powers between quasars and 
BL Lac objects, it can be concluded that most quasars have faster jets than most BL Lac 
objects. 

The second noteworthy inference that can be drawn from Figure El is that there is 
a continuous distribution of parsec-scale jet speeds going from BL Lacs to quasars. This 
supports the idea that the distinction between the BL Lac and quasar population, at an 
equivale nt width of 5A, is essenti ally an arbitrary one at least in terms of radio jet properties 



see also lScarpa fc Falomolll997l ) 



However, we note that inspite of the two blazar classes displaying a smooth continuation 
in properties (as is evident in Fig. El and other plots), they sometimes reveal differences 
in the correlation test results when considered separately In Table 12] we have listed the 
correlation test results separately for quasars and BL Lacs, if they differed from results for 
the combined blazar population. BL Lacs considered alone sometimes failed to exhibit the 
correlations observed in quasars (also see §3.1). Apart from the possible effects of small 
number statistics, a simple interpretation of this finding could be that unlike the quasars, 



-12- 



the BL Lacs do not c onstitute a homogeneous population. As has been previously s uggested 



in the literature (e.g.. lOwen et al.lll996l : ICara fc Lister! 120081 : lLandt fc Bignalll 120081 ). the BL 



Lacs might constitu t e the beamed population of both FRI and FRII radio galaxies. We note 
that iGiroletti et al.l (120041 ) have demonstrated that HBLs (which are absent in MOJAVE) 
conform to the standard unification scheme with respect to FRI radio galaxies. This again 
supports the proposed inhomogenity within the BL Lac class. 

Preliminary Monte Carlo simulations of a population of AGN similar to that of MOJAVE 
also show a correlation between apparent jet speed and extended luminosi ty (Cooper et al., 
i n pre paration). The simulations are based on the luminosity function of iPadovani fc Urrv 
(jl992j), and assume that the extended luminosity is unbeamed and proportional to the in- 
trinsic unbeamed parsec-scale luminosity. We are currently i ncorp orating the luminosity 
function derived from the MOJAVE sample by lCara fc Lister! ( 120081 ) into the simulations. 



3.4. Parsec-to-kiloparsec Jet Misalignment 



Apparent misalignm ent in the jet direction from parsec to kiloparsec scales is commonly 



observed in blazars (e.g.. |Pearson &: Readheadlll988l ). This misalignment could either be due 



to actual large bends in the jets, or due to small bends that are amplified by projection. 
We present the parsec- to kiloparsec-scale jet misalignment for the MOJAVE sources in 
Fi gure El The procedu re for estimating the parsec-scale jet position angles (PAs) is described 



in 



Lister et al.l (l2009al ). Kiloparsec-scale jet position angles were determined for the FRII 



sources using the brightest hot spots, especially when no jet was clearly visible, under the 
assumption that the brightest hot spot indicates the approaching jet direction due to Doppler 
boosting. The AIPS procedure TVMAXFIT was used to obtain accurate peak pixel positions 
of the core and hot spot, which were then used in the AIPS verb IMDIST to obtain the jet 
position angle. For the FRI-type sources we used TVMAXFIT for the core pixel position, 
and the AIPS verb CURVALUE to get the pixel position roughly towards the center and 
end of the broad jet/lobe. Finally IMDIST was used to obtain the jet position angle. Note 
that unlike regular radio galaxies, FRI blazars are typically one-sided, which makes the 
approaching jet direction easy to identify. Due to the wide jet/lobes in FRIs, the kiloparsec- 
scale jet position angle could be uncertain by 10° to 15°. When a jet feature appeared only 
partially resolved (viz., 1038+064, 1417+385), a two Gaussian component model was used 
in JMFIT to obtain pixel positions of the core and jet feature. 

We note that our approach (viz., of using the brightest hot spot to determine the 
approaching jet dir ection, when no jet was clearly visible) differs from that adopted by some 



other authors (e.g.. lXu et al.lll994l ) who assume that, when no jet is visible, the hot spot on 
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the side of the parsec-scale jet, represents the kiloparsec-scale jet direction. Our approach 
has led us to identify many more sources (~ 30%) with misalignment angles greater than 
90°. 




Fig. 6. — Misalignment angle between the parsec-scale and the kiloparsec-scale jet for the MOJAVE 
sources. There are signatures of a weak 90° bump. 



For an otherwise straight jet having a single bend, the apparent bending angle (77) is 
related to the intrinsic bending angle (ip) through the relation, cot 77 — { "' " M " '' 



where 9 is the angle to line of sight, and (f) is the azimuth of the bend ( lAppl et al.lll996l ). The 
fact that we see a majority of sources having close to zero apparent misalignment implies 
one of the following scenarios: (i) the intrinsic misalignment angle is small; (ii) both the 
intrinsic misalignment and jet inclination angles are large; (hi) the azimuth of the bend is 
close to 180°, i.e., the plane of the bend is perpendicular to the sky plane. Since scenario (ii) 
is unlikely for these blazars, possibilities (i) and (iii) are more favorable. In §4, we explain 
how the MOJAVE selection criteria could make it biased towards bent jets. 



Pearson fc Readheadl ( 119881 ) and lConway fc Murphyl ( 119931 ) reported an unexpected sec- 



ondary peak at APA = 90° in the misalign ment angle distribution. Th e signatures of a weak 
90° bump appear to be present in Figure [6j I Conway &: Murphyl ( 119931 ) concluded that while 
the largely aligned population could be explained by straight parsec-scale jets and small 
intrinsic bends between parsec- and kiloparsec-scales, the secondary peak sources must have 
gently curving (low-pitch) helical parsec-scale jets. This helical distortion could arise due 
to Kelvin-Hemlholtz instabilities or precession of the jet ejection axes. We reach a similar 
conclusion on the jet misalignments, following a different approach in §4.3. 
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DISCUSSION 



4.1. Selection Effects in MO J AVE 



The MOJAVE survey appears to select many sources with intermediate radio powers 
and radio morphology (e.g., high radio power BL Lacs, some with hot spots, and low radio 
power quasars). This could be a result of the MOJAVE selection criteria, which are based 
on the relativistically-boosted, high radio frequency (15 GHz) parsec-scale flux densities. 
This would make it more likely to encompass a much larger range in intrinsic radio powers, 
compared to other samples selected on the basis of total kiloparsec-scale flux-densities (e.g., 
th e 3CR sample) . The intrinsic radio luminosity function derived for the MOJAVE sample 
by ICara fc Listed (120081 ) also supports this view. 



Furthermore, as the MOJAVE survey picks up bright radio cores, most of the sources 
could either be fast jets pointed towards us, or curved/bent jets that have at l east some por- 
tion o f the jet aligned directly into our line of sight, at at least one epoch (e.g.. lAlberdi et al. 
19931 ). The MOJAVE survey could therefore be biased towards bent jets. Note that in the 
wide angle tail qua sars, the hot spots closer to the cores could be produced at the base 
of the plumes (e.g. , Hardcastle fc Sakellioul 120041 ). or indicate a sharp bend in the jet (e.g. 



Alberdi et al. 1993 



Jetha et al. 



200 



^). We discuss the possibility of beamed WAT quasars 



in the MOJAVE sample further in §4.3. 



4.2. Orientation Effects in the Sample 



Since the MOJAVE sample consists almost entirely of blazars, the sources are expected 
to have jets aligned close to our line of sight. Our examination of orientation effects in the 
sample, however, using two different statistical orientation indicators, R c and R v , reveals 
that a range of orientations must in fact be present in order to account for several observed 
trends. R c and R v also reveal, to some extent, dissimilar trends. We finally reach the 
conclusion that R v is a better indicator of orientation for this blazar sample. 

The ratio of the beamed radio core flux density (S core ) to the unbeamed extended radio 
flux density (S ext ), viz., the radio core prominence parameter (i? c ), ha s routinely been used 
as a statistical indicator of Doppler beaming and thereby orientation ( lOrr fc Brownelll982l ; 
Kapahi fc Saikialll982l : IrQiarb fc Shastril l2004f l . The fc-corrected R c (= f^f (1 + ^) CWe ~ aext , 
with a core = 0, a ext = 0.8) is plotted against the radio core luminosity in Figure UJ As 
expected, R c is correlated with the radio core luminosity (Table [2]). R c however, shows a 
significant anti-correlation with respect to the extended radio luminosity (Fig. [H]). R c does 
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Fig. 7. — The orientation indicators, R c (Left) and R v (Right) versus the 1.4 GHz core luminosity. 
Open and filled circles denote quasars and BL Lacs, respectively, while open squares denote radio 
galaxies. Core-only sources are represented as upper limits. 



not show a correlation with jet misalignment (Table [2]). In fact the quasars considered alone 
showed a weak anti-correlation between R c and jet misalignment (implying that the more 
core-dominant sources have smaller jet misalignments). These results are unexpected, be- 
cause in the former case, the extended radio luminosity is expected to be largely unbeamed, 
while in the latter, the jet misalignment is affected by projection and therefore orienta- 
tion (in the sense that the more cor e-dominant source s show larger jet misalignments, e.g., 
Kapahi fc Saikialll982l ). As an aside. iBrothertonl (119961 ) have noted that the anti-correlation 
between logi? c — \ogL ext may be suggesting that core-dominated quasars are intrinsically 
fainter than their lobe-dominated counterparts. This would be consistent with the idea that 
the MOJAVE sources span a large range in intrinsic radio powers (cf. §4.1). 

We also examined the relationship between R c and apparent jet speed. While both 
R c and f3 app depend on the Lorentz factor and orientation, many of the viewing angles 
may be inside the critical angle for maximum superluminal speed, which would spoil any 
correlation. However, the quasars considered alone do show a significant anti-correlation, 
contrary to expectations. But as we see below, the alternate orientation indicator, R v , does 
show the expected behavior. 



Wills fc Brothertonl (119951 ) defined R v as the ratio of the radio core luminosity to the k- 
corrected absolute V-band magnitude (M a b s ): logR v = logh^- = (logL core +M a b s /2.5) — 13.7, 



where M a b s = My — k, and the /c-correction is, k = —2.5 log (1 + z) 1 a ° pt with the optical 
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Fig. 8. — The orientation indicators, R c (Left) and R v (Right) versus the 1.4 GHz extended 
luminosity. Open and filled circles denote quasars and BL Lacs, respectively, while open squares 
denote radio galaxies. Core-only sources are represented as upper limits. 



spectral index, a opt = 0.5. R v is suggested to be a better orientation indicator then R c 
since the optical luminosity is likely to b e a better measure of intrinsic jet power (e.g., 
Maraschi et al.l 120081 ; iGhisellini et al.l 120091 ) than extended radio luminosity. This is due to 
the fact that the optical contin uum luminosity is correlated with the emission-line luminosity 
over four orders of magnitude (jYee fc Okdll978T). and the emission-line luminosity is tightly 
correlated with the total jet kinetic power feawlings fc SaunderslllQQll L The extended radio 
luminosity, on the other hand, is suggested to be affected by interaction with the environment 
on kiloparsec-scales. 

Figure [7] suggests that R v is indeed a better indicator of orientation as the correlation 
with radio core luminosity gains in prominence (see Table |2j). R v shows a weak positive 
correlation with the extended radio luminosity (Fig. [8]). Partial regression analysis however, 
shows that the linear correlation between logi?„ and logL ext is no longer significant when the 
effects of \ogL core are removed (r X y.z = 0.078, t statistic = 0.84, p = 0.4013). This implies 
that the extended radio emission is largely unbeamed. R v is correlated with the parsec-to- 
kiloparsec jet misalignment, suggesting that orientation does play a role in the observed jet 
misalignments, at least for the quasars. The lack of such a correlation for the BL Lacs could 
again be interpreted to be consequence of their comprising of not only beamed FRIs (which 
could have intrinsically distorted jets) but beamed FRIIs as well. We discuss this point again 
in §4.5. Finally, there appears to be no correlation between R v and (3 a pp, as is expected if 
many blazar viewing angles are inside the critical angle for maximum superluminal speeds, 
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which would then result in smaller apparent speeds. This result is again consistent with R v 
being a better indicator of orientation than R c . 

Here we would like to note that the correlation between apparent jet speed and absolute 
optical magnitude does not seem to be as significant as the one observed between apparent 
jet speed and extended radio luminosity (Table 2). This could in principle undermine the 
suggestion that the optical luminosity is a better indicator of jet kinetic power than extended 
radio luminosity However, since the optical luminosity is more likely to be affected by strong 
variability, than the extended radio luminosity, the lack of a strong correlation can perhaps be 
accounted for, in these non-contemporaneous radio and optical observations. Nevertheless, 
this is an important result to bear in mind, that requires further testing. 



4.3. Environmental Effects 
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Fig. 9. — (Left) 1.4 GHz extended radio luminosity versus the absolute optical magnitude. (Right) 
The environment indicator versus redshift. Open and filled circles denote quasars and BL Lacs, 
respectively, while open squares denote radio galaxies. Core-only sources are represented as upper 
limits. 



It has been suggested that apart from projection and relativistic beaming effects, the 
complex radio structures in quasars could arise due to inte ractions with the surrounding 
medium (IPearson fc Readheadl Il988l ; iBarthel &: Mileyl Il988l ) . As is observed in Figure [91 
the extended radio luminosity appears to be correlated with absolute optical luminosity. 
However, a partial regression analysis shows that the linear correlation becomes weak or non- 
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existent, after the effects of luminosity distance are removed (for quasars, r X y.z = —0.169, 
t statistic = -1.62, p = 0.1088; for BL Lacs, r XY . z = 0.079, t statistic = 0.31, p = 0.7608). 
The lack of a strong correlation may again be attributable to strong optical variability 
in these blazars. If the extended radio luminosity is indeed affected by interaction with the 
kiloparsec-scale environment, and the optical luminosity is closely related to the AGN power, 
the ratio [log L PX t ]/M nh , q , can serve as a pro be for environmental effects on kiloparsec-scales 
as suEre; ested by Iwills fc BrothertorJ Il995h . Henceforth we will refer to this ratio as the 



"environment indicator" (we will drop the "log" term in the ratio for convenience, in the 
following text). 

A higher value of the ratio could indicate greater jet-medium interaction, (as say, from 
an asymetrically dense source environment). Dense environments can decrease expansion 
losses in the source, but increase radiative losses, making the sources brighter at low radio 
frequencies (e.g., in Cygnus A, see iBarthel fc Arnaudlll996h . Note that the optical magnitude 



in blazars is expected to be related to the AGN (accretion disk or jet), rather than to 
starlight. This might not be true for the nearby radio galaxies, where starlight might indeed 
be a significant contributor to the optical magnitude. The L ext /M a t> s ratio might therefore 
not be a good environment proxy for the radio galaxies. We have therefore excluded radio 
galaxies from the correlation tests related to the environment proxy indicator. 

We plot the environment proxy ratio with respect to redshift in the right hand panel 
of Figure M and find a strong positive correlation. The "alignment effect" between the radio 
source axis and the emission line region in high redshift radio sources demonstrat es that local 



galactic asymme tries (on scales of several kiloparsecs) increase with redshift (lMcCarthvlll993l : 



Best et al.lll999l ). This result therefore, is consistent with the suggestion of the environmental 
dependence of the extended luminosity. We note that the BL Lacs considered alone fail to 
show a strong correlation (Table which could be a consequence of an inhomogeneous 
parent population, as discussed in §3.3. 

We find that while the parsec-scale jet speeds appear to be tightly correlated with 
the extended radio luminosity (Fig. E]), and weakly correlated with the optical luminosity, 
they seem to not show any correlation with the kpc-scale environment indicator (Fig. [TO]) . 
This could imply that the parsec-scale jet speeds are not driven by the kiloparsec-scale 
environment, but rather by factors intrinsic to the AGN, or the AGN's local parse c-scale 



envir onment. For instance, the jet speeds could be related to the black hole spins (e.g., [Meier 



19991 ) . However, like other results pertaining to optical luminosity, this lack of a correlation 



needs further testing with contemporaneous optical and radio data. 



We find that the parsec-to-kiloparsec jet misalignment is not correlated with redshift 
(Fig. ITT]) . This goes against the idea that the bending of the jets is influenced by interaction 
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Fig. 10. — Apparent parsec-scale jet speed versus the environment indicator. Open and filled 
circles denote quasars and BL Lacs, respectively, while open squares denote radio galaxies. 



with the medium; an effect which should g et more prominent w ith i ncreasing re d shift. This 
result is c onsistent with t he su ggestions of Hintzen et al. ( 1983 ) and Appl et all fll996l ). We 
note that iHintzen et al.l ( 19831) used the kiloparsec-scale jet-to-counterjet misalignment in 
their study. iKharb et al.l (l2008bl ) found a similar lack of correlation between jet-to-counterjet 
misalignment and redshift for an FRII radio galaxy sample. Thus, the jet misalignment could 
also be related to factors closely associated with the AGN itself. For instance, the presence of 
binary black holes or kicks imparted to the bl ack hole via a black hole merger, could switc h 
the ejection direction on a short timescale (e.g.. lBegelman et al.lll980l ; iMerritt fc Ekersll2002l ). 
Black hole re- alignment due to a warped accretion disk could also give rise to jet precession 
(lPringldll997l ). We note that a few sources have morph ologies resembling fast, precessing jets 
("banana" shaped) as modelled by iGower et al.l (119821 ). Alternatively, the jet misalignment 
could be influenced by the AGN's local parsec-scale environment. 

The jet misalignment is inversely correlated with the kpc-scale environment indicator 
(see right hand panel of Fig. ITTi) . The inverse relation suggests that jets with smaller mis- 
alignments (i.e., straighter jets) show signatures of greater environmental effects on their lobe 
emission. This could either be suggestive of a uniformly dense confining medium on different 
spatial scales around the source, or, of projection effects in the sample, which could boost the 
optical emission and decrease the L ext /M a t> s ratio in sources with larger jet misalignments. 
If the latter is true, then this is an important caveat — the environment indicator could 
also be affected by source orientation in blazars. However, the orientation-dependence alone 
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Fig. 11. — (Left) Jet misalignment versus the environment indicator. (Right) Misalignment angle 
versus redshift. Open and filled circles denote quasars and BL Lacs, respectively, while open squares 
denote radio galaxies. 



cannot satisfactorily explain the correlation between the environment indicator and redshift 

(Fig. ED. 

Many MOJAVE blazars exhibit radio structures resembling those of wide angle tail 
radio galaxies. Obs ervations have indicated that WAT quasars could be mem bers of rich 
galaxy clusters (e.g., Hintzen et al. 1983 : Harris et al. 1983 : Blanton et al.ll200ol ). Using the 
NASA/IPAC extragalatic database (NED) we found that nearly 7% of the MOJAVE sources 
have a galaxy cluster ~ 5' away, while 14% have a galaxy cluster ~ 10' — 15' away. While 
this fraction appears to be small, we must keep in mind that most of the blazars are at high 
redshifts, while the cluster information is highly incomplete even at redshifts z > 0.1 (e.g., 



Ebeling et al. 



1998 



Bohringer et al.ll2004l ). Thus the possibility remains that a large number 



of MOJAVE blazars inhabit galaxy clusters. This suggestion has interesting ramifications, 
as discussed below. Furthermore, this could help us understand some of the results discussed 
in this section. 



4.4. Relation Between Radio Power and Emission-line Luminosity 



Rawlings fe Saunders! (Il99ll ) demonstrated a tight correlation between emission-line lu- 
minosity and total jet kinetic power in a sample of radio-loud AGNs. This correlation 
was tighter than the one observed between emission-line luminosity and radio luminosity. 



- 21 - 



Landt et al.l (120061 ) discovered a large number of "blue quasars" with strong emission lines, 
but relatively low radio powers in the Deep X-ray Radio Blazar Survey (DXRBS) and the 
RGB samples, and suggested that this went against the close relationsh ip between emission- 
line and jet power. These and other observations (e.g., Ku et al.lll999l ) could be suggesting 
that the mechanism for the production of the optical- UV photons that ionise the emission-line 
clouds is decoupled from the mechanism that produces powerful radio jets. Many MOJAVE 
blazars also seem to not follow the relation between emission-line and radio luminosity; that 
is, some quasars and BL Lacs have relatively low and high radio powers, respectively, and the 
suggestion of decoupled mechanisms might hold true. However, an alternate hypothesis by 
which their emission-line luminosity could still be correlated with their intrinsic jet kinetic 
power, could come about if the outlying BL Lac objects were present in dense confining en- 
vironments, while the outlying quasars were old, field sources. Confinement would increase 
the synchrotron radiative loss es due to the amplificatio n of the magnetic field strength, mak- 
ing the sources brighter (see iBarthel fc Arnaudlll996l ). while expansion losses and gradual 
radiative losses could presumably reduce the low frequency radio emission in the old, field 
quasars. 

It has been suggested t hat a large fractio n (~50%) of the high redshift radio sources 



could be GPS sources (e.g., lO'Dea et al.lll99ll ). which are likely to be present in confined 
environments. One of the highest redshift MOJAVE sources, 0742+103 (z = 2.624), is a 
GPS quasar. There are several other MOJAVE sources that have peaked spectra, but do not 
show the low-variability or the double-lobed parsec-scale structures that are characteristic 
of the GPS class. These could be "masquerading GPS" sources, whose overall spectrum 
happens to be dominated by an unusually bright feature in the jet, located downstream 
from the core. 



4.5. Blazar Division Based on Radio Power 

In order to gain more insight into the primary factors that drive the FR dichotomy, 
in this section we disregard the emission-line division (i.e., quasars, BL Lacs), and divide 
sources into FRI, FRI/II and FRII classes based solely on their 1.4 GHz extended powers 
(following Fig. HJ). Figure [T2l plots the environment indicator, L ext /M a f, s , with respect to 
redshift, using the new classification. A quick comparison with Fig. M reveals that many 
of the low redshift quasars fall into the intermediate FRI/II class. The low redshifts of 
these sources discounts the possibility of a (1 + z) 4 surface brightness dimming effect in 
them, which could potentially reduce their extended luminosity, making them fall in the 
intermediate luminosity class. This finding is consistent with the discovery of "blue quasars" 
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by lLandt et al.l (120061 ). which have radio powers and spectral energy distributions (SEDs) 
similar to the high energy peaked BL Lacs. 
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Fig. 12. — The environment indicator versus redshift, with sources divided on the basis of extended 
radio luminosity into FRI (=1), FRII (=2), and FRI/II (=filled circles). 

Most beamed FRIs lie at redshifts below 0.5 in Fig. [12j This could be a limitation 
of the flux limited sample. Interestingly, the FRI/II sources lie at all redshift ranges, and 
throughout the L ext /M abs distribution. This suggests that these "intermediate" luminosity 
sources are not produced as a result of different environmental conditions compared to, say, 
other FRIIs. Similarly, sources divided on the basis of the presence or absence of hot spots, 
show no discernable difference in the L ex t/M a b s vs. z pl ane. It is important to note that 
environmental radio boosting, as mentioned in §4.4 (also iBarthel fc Arnaudl 119961 ) could be 
instrumental in blurring the Fanaroff-Riley dividing line in the MOJAVE blazars. 

Lastly, the orientation indicator, R v , seems to be correlated with misalignment only 
for FRIIs, but not for FRIs or FRI/IIs (Table [2]). This implies that projection is playing a 
significant role in the jet misalignment of FRII sources, but the jets are intrinsically more 
distorted in FRIs and FRI/IIs. A corollary to this is that the differences between FRIs 
and FRIIs could be related to factors affecting the black hole spin direction, since we have 
previously concluded that jet mislignments appear to be largely unrelated to kiloparsec-scale 
environmental effects. 
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4.6. Extremely Misaligned Sources 



Six MOJAVE sources (viz., 0224+671, 0814+425, 1219+044, 1510-089, 1828+487, 
2145+067) have jet misalignments greater than 160°. Such extremely misaligned sources 
could be bent jets viewed ne arly head-on, as suggested for the gamma-ray blazar PKS 
1510—089 ( iHoman et al.ll2002l ). An alternate sc enario could be that some of th ese misaligned 
sources are "hybrid" morphology sources (e.g., lOopal-Krishna fc Wiitall2~000n . Three of six 
sources with large misalignments (viz., 0224+671, 1219+044, 1510—089) indeed seem to pos- 
sess radio morphologies that could be classified as "hybrid". The large misalignment then 
makes sense: the VLBI jet (used to determine the parsec-scale jet position angle) could be 
on the FRI side, while the sole hot spot (used to determine the kiloparsec-scale jet position 
angle) could be on the FRII side. 

One interesting source that shows extreme misalignment (~ 135°) is M87 (1228+126) . 



M87 appears to e xhibit a weak hot spot on the counterjet side (e.g., lOwen et al.l Il990l ; 
Sparks et al.lll992l ). Since the brightest emission on the counterjet side is at the extremity 
of the radio lobe (like in FRIIs), M87 could be classified as a hybrid radio morphology 
source. I ts radio luminosity ( L 17s = 1 x 10 2 5 W m ~ 2 ) also places it close to the FRI/FRII 
division ( Biretta et al.l fl999l ) . iKovalev et al. ( 2007 ) have observed a very short, weak radio 
counterjet on parsec-scales in M87, but perhaps that is just a slower outer sheath that quickly 
terminates. 

If indeed some of the extremely misaligned sources are "hybrid" sources, the lack of a cor- 
relation between jet misalignment and environment appears to contradict t he suggestion tha t 
asymmetries in galactic environments, whic h increase with redshift (e.g., iBest et al.111 999) 



give rise to hybrid morphology sources (e.g., iGopal-Krishna fc Wiitall2000l ; Miller fc Brandt 



20091 ) . Moreover, the fraction of such sources would be higher in the M OJAVE yoyo — 070J 
than that reported for the FIRST survey (1%, iGawrohski et al.l 120061 ) . suggesting that hy- 
brid sources may not be as rare as previously supposed. It is also interesting to note that 
not many MOJAVE "hybrid" morphology sources fall in the FRI/II luminosity range, but 
rather have FRII luminosities. In other words, the MOJAVE "hybrid" morphology sources 
are not necessarily "hybrid" in terms of radio power, but mostly fall into the FRII luminos- 
ity class. Chandra X-ray observations of these sources could resolve some of these apparent 
contradictions. 
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4.7. Core-only Sources 



About 7% of the MOJAVE sources do not show any extended radio emission in our 
1.4 GHz VLA A-array images. However the ratio of their 5 GHz radio flux density to their 
B-band optical flux density is of the order of a few hundred, placing them firmly in the 
radio-loud AGN class 



(e.g. 



Kellermann et al.l Il989l ) . The dynamic range achieved in the 



radio images of these sources varied from about 6000:1 (e.g., 0955+476) to 30,000:1 (e.g., 
0202+149), being typically of the order of 10,000:1. While some sources, like 1749+096, failed 
to re veal any extended emission in a number of different datasets (see also iRector fc Stocke 
20011 ). we did detect faint hal o -like e mission around 1548+056, which was listed as a core- 
only source by Murphy et al.l ( 119931 ). Note that observations with the VLA A-array-only 
are likely to miss very diffuse radio structure, due to the lack of short baselines. We believe 
therefore that these core-only sources probably have associated faint emission which require 
more sensitive observations to detect. 

The observed fraction of c ore-only sources in MOJAVE is smaller than t hat reported 
for other quasar surveys (20%, iGower fc Hutchingslll984J : IPadrielli et al.l 119881 1 . While the 
core-only sources could have been affected by the redshift surface brightness dimming effect, 
our discussion in the previous section argues against this conjecture. Alternatively, these 
sources could just be normal quasars at very small angles to line of sight, so that the VLA 
resolution is insufficient to delineate the various components. This appears to be the case in 
1038+064, and 1417+385, where a two Gaussian component model fits the "core" emission 
better than a single one. The hypot hesis that at least som e of these quasars could be beamed 
FRIs also cannot be ruled out (e.g.. iHey wood et al.l 120071 ) . 



4.8. Optical Polarization Subclasses 



Low optical polarization radio quasars (LPRQs) consistently reveal core optical frac- 
tional polarization, m opt < 3%, while the high optical polarization quasars (HP Qs) routinely 
reveal highly polarized optical cores with m opt > 3% (jAngel fc Stockman! Il980h . The classi- 
ficati on of the MO J AVE sourc es into these categories is listed in Tabled! and were obtained 
from llmpey et al.l (1199 ll ) and iLister fc Smith! (120001 ). In all, there are 38 HPQs and 16 
LPRQs in the MOJAVE sample. We find no statistically sign ificant differences in the core 
or extended radio luminosity between LPRQs and HPQs. In iKharb et al.l ( 12008al ) we had 
noted a difference in the integrated (single-dish) 1.4 GHz radio luminosity between LPRQs 
and HPQs for redshifts less than one. We do not find a difference in the total radio lumi- 
nosity using the present higher resolution VLA observations. Furthermore, we do not find 
a statistically significant difference in jet misalignment between LPRQs and HPQs. The 
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parsec-scale apparent jet speeds of HPQs, however, are systematically higher than those 
in LPRQs (Kolmogorov-Smirnov test probability that the speeds in HPQs and LPRQs are 
similar, is p=0.03). Additional optical polarimetry information on the sample is needed to 
further investigate these trends. 

5. SUMMARY AND CONCLUSIONS 

1. We describe the results from a study of the extended emission at 1.4 GHz in the 
MOJAVE sample of 135 blazars. New VLA A-array images of six MOJAVE quasars, 
and previously unpublished VLA A-array archival images of 21 blazars, are presented. 
The kiloparsec-scale structures are varied, with many sources displaying peculiar, dis- 
torted morphologies. Many blazar jets show parsec-to-kiloparsec-scale jet misalign- 
ments greater than 90°. These characteristics have been reported in other quasar sur- 
veys as well. The 90° bump in the jet misalignment distribution (of which we observe 
a weak signature) has been suggested to be a result of low-pitch helical parsec-scale 
jets in the literature. 

2. A substantial number of MOJAVE quasars (~ 22%) and BL Lacs (~ 23%) display radio 
powers that are intermediate between FRIs and FRIIs. Many BL Lacs have extended 
luminosities (~ 27%) and hot spots (~ 60%) like quasars. The hypothesis that at least 
some of the core-only quasars are in fact FRI quasars, cannot be discarded. In terms 
of radio properties alone, it is difficult to draw a sharp dividing line between BL Lacs 
and quasars, in the MOJAVE sample. This could be a result of the MOJAVE selection 
effects which naturally pick sources with a large range in intrinsic radio power (low 
power quasars and high power BL Lacs), and preferentially bent jets. While the quasars 
and BL Lacs display a smooth continuation in all their properties, the correlation test 
results differ sometimes when they are considered separately. This can be understood 
if the BL Lacs did not constitute a homogeneous population like the quasars, but had 
both FRI and FRII radio galaxies as their parent population. These findings challenge 
the simple radio-loud unified scheme, which links FRI sources to BL Lac objects and 
FRIIs to quasars. 

3. We find that there is a significant correlation between extended radio luminosity and 
apparent parsec-scale jet speeds. This implies that most quasars have faster jets than 
most BL Lac objects. The large overlap between many properties of quasars and 
BL Lacs (e.g., 1.4 GHz radio power, parsec-scale jet speeds) however suggests that, at 
least in terms of radio jet properties, the distinction between these two AGN classes, 
at an emission-line equivalent width of 5A, is essentially an arbitrary one. 
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4. These observations suggest that the mechanism for the production of optical-UV pho- 
tons that ionise the emissi on-line clouds i s decoupled from th e mechanism that produces 
powerful radio jets (e.g., Xu et al. 1999 : Landt et al.lhood ). An alternate hypothesis 



could be that, BL Lacs with high radio powers are present in dense confining envi- 
ronments, which would increase the radiative losses in them and make them brighter, 
while the radio quasars with low radio powers are old, field sources. Environmental 
radio boosting ( Barthel fc Arnaudlll996l ) could also explain why the Fanaroff- Riley di- 



viding line is blurred in the MOJAVE sample. X-ray cluster studies could be used to 
test these ideas. 

5. The ratio of the radio core luminosity to the ^-corrected optical luminosity (R v ) ap- 
pears to be a better indicator of jet orientation than the traditionally used radio core 
prominence parameter (R c ). This seems to be a consequence of the environmental 
contribution to the extended radio luminosity, used in R c . Trends with R v reveal that 
even though the sample comprises largely of blazars, there is a reasonable range of 
orientations present in the sample. 

6. Trends with the "environment" proxy indicator, L ext /M a t> s , reveal that jet speeds seem 
to not depend on the environment on kiloparsec scales, while the jet misalignments seem 
to be inversely correlated. The jet misalignments are also not correlated with redshift. 
It appears that the parsec-to-kiloparsec jet misalignment, the parsec-scale jet speeds 
and to an extent the extended emission (which is related to jet speed) are controlled 
by factors intrinsic to the AGN. Black hole spins could dictate the jet speeds, while the 
presence of binary black holes, or kicks imparted to black holes via black hole mergers 
could influence the jet d irection. Thi s is con sistent with radio morphologies similar to 



precessing jet models of iGower et al.l (119821 ) present in the MOJAVE sample, and the 



signature of the 90° bump in the jet misalignment distribution, which is attributed to 
low-pitch helical parsec-scale jets. 

7. If some of the highly misaligned sources (jet misalignment ~ 180°) are "hybrid" FRI+II 
morphology sources, then the fraction of such sources is higher in the MOJAVE surve y 
Vo — 8%), than that reported for the FIRST survey (1%, iGawrohski et al.l 120061 ) . 



Furthermore, the lack of a correlation between jet misalignment and the environmental 
indicator would appear to contradict the suggestion that hybrid morphology is a result 
of the jet interacting with an asymmetric environment. X-ray observations of hybrid 
sources are needed to resolve these issues. 

8. About 7% of the MOJAVE quasars show no extended radio emission. We expect 
more sensitive radio observations to detect faint emission in these sources, as we have 
detected in the case of 1548+056, previously classified as a core-only quasar. The 
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hypothesis that at least some of these quasars could be beamed FRIs cannot be ruled 
out. 

9. While the extended radio power or jet misalignments do not show statistically signif- 
icant differences between the two quasar optical polarization subclasses, viz., LPRQs 
and HPQs, the parsec-scale apparent jet speeds of HPQs are systematically higher 
than those in LPRQs. 
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Table 1. The MO J AVE sample 
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Note. — Col. 1 & 2: IAU B1950 names of the MOJAVE sources and their redshifts r — sources rln not .shnm any discernable extended 

emission. Col. 3: Apparent V-band magnitude obtained from the AGN catalog of IVeron-Cettv fc Veronl (|2006h or NED. Col. 4: Blazar 
classification — Q = quasar, B = BL Lac object, G = radio galaxy, U = Unidentified. Col. 5: Highly optically polarised — Y T|^ PS '|N = 
No, .... = no data available. Col. 6 &; 7: Peak and rms intensity in radio map, respectively. I r ms for the sources in[ Murp hy et al.l (1993) are 
the CLEV values in their Table 2. J for Cygnus A, the peak intensity corresponds to the south-eastern hot spot. Col. 8: Central observing 
frequency in GHz. Col. 9: Apparent parsec-scale speed of the fastest jet feature. Col. 10: Integrated core flux density at ~1.4 GHz obtained 
with JMFIT. Col. 11: Extended flux density obtained by subtracting the core from the total radio flux density. Col. 12: Radio morphology — 
c = core only, ch = core-halo, lhs = 1 hot spot, 2hs = 2 hot spots, 1 = 1 sided, 2 = 2 sided, lhs+2 = 2 sided structure with hot spot on one 
side. Alternate morphologies are listed separated by commas. Col. 13: Parsec-scale jet position angle. Col. 14: Kiloparsec-scale jet position 

angle CxJ 1 ^-References and VT A amhive Pr^jg'-t IDs for the radio data — f=Ncw data presented in this paper. *=Data presented in 

ICooper et al.l i20071) . Murphy = lMurphv et al.l (U993T) . 
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Table 2. Correlation Results 
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Note. — Cols. 1 & 2: Parameters being examined for correlations. Results from a partial regression 
analysis are cited in the main text. The ten core-only sources have been excluded from correlations 
with extended luminosity. Including the upper limits in extended luminosity as detections, does not 
alter any of the observed trends. Col. 3: The correlation test results have been presented separately 
if they were different for quasars and BL Lacs considered alone. ALL = quasars + BL Lacs + radio 
galaxies, Q = only quasars, B = only BL Lacs, Q+B = quasars + BL Lacs. FRII, FRI/II, and FRI arc 
as denned in §4.5. Cols. 4 & 5: Spearman rank correlation coefficient and chance probability. Cols. 6 
& 7: Kendall Tau correlation coefficient and chance probability. Col. 8: Indicates if the parameters are 
significantly correlated. YES?' indicates a marginal correlation. We note that at the 95% confidence 
level, 1.7 spurious correlations could arise from the Ri35 correlations that we tested. 
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Fig. 13.— New radio images from the VLA archive of (a) 0727-115, (b) 0736+017, (c) 0742+103, 
(d) 0838+133, (e) 1124—186, and (f) 1334—127. The contours are in percentage of the peak surface 
brightness and increase in steps of 2. The lowest contour levels and peak surface brightness are, 
(a) ±0.01, 3.15 Jy beam" 1 ; (b) ±0.042, 2.33 Jy beam^ 1 ; (c) ±0.021, 3.61 Jy beam" 1 ; (d) ±0.17, 
289.3 mJy beam" 1 ; (e) ±0.085, 651.2 mJy beam" 1 ; (f) ±0.042, 1.94 Jy beam -1 . 
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Fig. 14 — New radio images from the VLA archive of (a) 1504-166, (b) 1548+056, (c) 1726+455, 
(d) 1733-130, (e) 1741-038, and (f) 1751+288. The contours are in percentage of the peak surface 
brightness and increase in steps of 2. The lowest contour levels and peak surface brightness are, (a) 
±0.042, 2.36 Jy beam" 1 ; (b) ±0.021, 2.09 Jy beam" 1 ; (c) ±0.042, 998.9 mJy beam" 1 ; (d) ±0.021, 
6.11 Jy beam" 1 ; (e) ±0.042, 1.68 Jy beam" 1 ; (f) ±0.17, 263.5 mJy beam" 1 . 
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Fig. 15.— New radio images from the VLA archive of (a) 1807+698, (b) 1936-155, (c) 1958-179, 
(d) 2005+403, (e) 2008-159, and (f) 2021+317. The contours are in percentage of the peak surface 
brightness and increase in steps of 2. The lowest contour levels and peak surface brightness are, 
(a) ±0.042, 1.14 Jy beam" 1 ; (b) ±0.085, 1.05 Jy beam" 1 ; (c) ±0.085, 1.77 Jy beam" 1 ; (d) ±0.085, 
2.43 Jy beam" 1 ; (e) ±0.085, 544.5 mJy beam" 1 ; (f) ±0.085, 2.97 Jy beam" 1 . 
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Fig. 16. — New radio images from the VLA archive of (a) 2021+614, (b) 2037+511, and (c) 
2128—123. The contours are in percentage of the peak surface brightness and increase in steps of 2. 
The lowest contour levels and peak surface brightness are, (a) ±0.042, 2.64 Jy beam -1 ; (b) ±0.042, 
4.61 Jy beam" 1 ; (c) ±0.042, 1.37 Jy beam" 1 . 



